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Q^ ' Abstract 

Consistent quantum formalism based on the localized basis of the Wannirer functions in Heisen- 
O ' berg and Schrodinger pictures to describe propagation of electromagnetic field in a three dimen- 

C^ [ sional media including diffraction is presented. In the Schrodinger picture the Fokker-Planck equa- 

tion for the Glauber-Sudarshan quasiprobability and corresponding Langevin equations are given. 
As result the space-time description is obtained by a simple changing variables in the temporal 
master equation of the field. Using this formalism it is shown that the existence of integrals of 
^ I motion in the propagation of light in a medium under the condition of nondegenerated parametric 

CO ■ and two-photon interactions results in amplification of modes when nonclassical properties of the 

light are conserved. Quantum propagation of light in a linear medium taking into account the 
diffraction is considered and its solution is found. 
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1 Introduction 

When considering statistics of the light the usual approaches are based on the temporal evolution of 
the electromagnetic field. It is natural for the problems when the spatial behavior can be neglected as 
for the high-Q optical cavity scheme inside which one the temporal feature is important. By contrast 
the spatial behavior plays the key role in propagation of the light through media or distributed system 
that can be not considered as a point one. In these systems new states of the light named the spatially 
squeezed states arise O. The quantum formalism for distributed systems has been developed in a 
number of works. One-dimensional approaches were presented in Q. Theory for continuous variables 
in the Heisenberg picture was given in [^, where the continual coherent states and squeezed states of 
the field have been introduced and theory of the light detection in free space has been presented. With 
the use of these approaches evolution of the light statistics has been examined for the three-photon 
Q and four-photon ^ parametric interaction in transparent media. The multi-photon processes was 
discussed in Q. 

In this paper we present a quantum formalism which enables to consider propagating of the light 
in a three-dimensional medium with diffraction to be taken into account. It is based on a localized 
basis of the Wannirer functions, and allows the conversion of the temporal description to space- 
time description by a simple change of variables. Indeed it is the way used in solid state physics 
for transition from collective to local variables Q . The main feature of the local description in the 
Wannirer basis is that the problem turns out to be multi-particle when the local field oscillators 
interact among themselves even in free space. From the physical point of view this is a transfer of 
excitation or the light propagation process. In other hand the interaction of local oscillators produces 
the coupled equations which as we show have the form of the transfer equations. Our approach is 
formulated in both the Heisenberg and Schrodinger pictures where the Fokker-Planck equations for 
the field quasiprobability are derived together with the corresponding Langevin equations. 

To describe interaction of light with atoms the adiabatic illumination of fast atomic variables is 
often used to obtained a closed equation of field or master equation |§], § that is the start point 
to analyze statistics of the light. The simple recipe of how to get the local description allows to 
obtained the transfer equation immediately from the temporal master equation missing all steeps of 
the derivation procedure. 

The potential of the presented formalism is illustrated by two problems. In the first one the 
statistics of the light is considered for the propagation in a one-dimensional medium with two-photon 
and parametric interaction. Here within the framework of local description an interesting peculiarity 
due to the integrals of motion arises. The existence of the integrals makes it possible to establish 
the main features of the statistics transformation immediately without going into the solution of the 



dynamic problem [10|. It has been shown by this way that enhancement of light and the conservation 
of nonclassical properties can be possible. In the second problem the propagation of the light in a 
linear three-dimensional medium with diffraction is considered and the solutions have been found. 

2 Wannirer basis. Local operators 

In one-dimensional normalization space L plane waves ipk{x) = (1/vL) exp(zA:x), where A/c = 2it/L, 
are forming a complet orthonormal basis. It is used for a standard representation of the electromagnetic 
field strength operator, where the operators of photon creation and annihilation al, a^ are arising with 



commutational relations 

[ak,al,] = 6kk' (2.1) 

Operators al,,ak describe the creation and annihilation of photons of a wave vector k in the whole 
space L. These operators may be named collective, because they are responsible for the excitation of 
the whole volume. For the local description of the electromagnetic field instead of plane waves we use 
Wannirer functions known in solid state physics Q , which are packets of plane waves 

Wm{x -l) = —j= V ex.Y>{-ikl)ipk{x) (2.2) 

V-ZV , 

the wave vectors k ^ m are lying in a band m — vr/a < k < m + tx ja. Here each band m, or 
zone, corresponds to the partition of one-dimensional space L into N = L/a cells, which positions 
are determined as / = no. In k-space the centers of so defined bands are separated by an interval 
Am = 27r/a. The terms in ( |2.2| ) are forming a geometrical progression, so summation is possible. Then 
it easily verified that for large N the Wannirer functions are localized in a cell with a coordinate I in 



region ~ a. The Wannirer functions defined in accordance with (2^) form the complet orthonormal 
set 

^Wm{x-l)w*^{x' -I) = 5{x-x') (2.3) 

ml 

{Wm{x -l),W*^,{x-l')) = 6mm'h' (2.4) 

where the scalar product is defined by the integration over the whole space L. 

Wannirer basis serves as a framework for the introduction of local operators of creation and anni- 
hilation a}^i,am{l): 

^ipk{x)ak = ^am{l)wm{x -I) (2.5) 

k ml 

From ( p. 5]) it follows that the operators are coupled by a unitary transformation 

a^{l) = Y.C*^,{l)ak (2.6) 

k 

Ofc = X] Cmk{l)am{l) (2.7) 

ml 

where 

Cmk{l) = ^ E ^M-ik'l)5k'k (2.8) 

V ^^ k'^m 

For local operators the following commutational relations are valid: 

[0'm{l),a)rn'^^')\=^rnm'^ll' (2.9) 



The relations ( |2.5| ) allow to interpret them as operators of creation and annihilation of a photon in 
a point / in the vicinity ~ a. As a result local field oscillators described by the introduced operators 
are defined in the space L. Prom here on we will refer the packets of plane waves forming the local 
operators as local modes with wave number m and with a width Az/^ = c/a. Strictly speaking these 



packets are not modes, because there are no physical reasons for their distinctions, however the term 
seems to be convenient. 

The equation ( |2.5D is a starting point in the formahsm of quantum transfer theory, which may be 
formulated in different pictures. From the operational point of view the unitary transformations ( ^.6| ) 
and ( |2.7D following from ( p.5| ) mean that the transition from the nonlocal to local description and vice 
versa is accomplished by the change of variables. 

The specific feature of local description is its many-particle character, where the local oscillators 
are interacting already in the free space. This interaction describes the excitation transfer or light 
propagation process, which determine the structure of equations of motion, where chains of the BBGKI 
type (Bogolyubov, Born, Green, Kirkwood, Ivon) are arising resulting finally in the propagation 
equations. 

Consider the case of free space, for which the field evolution is defined only by a Hamiltonian 
Ho = J2k 'hoJka\ak. Uk = ck. Let the operator a^ in (|2.5|) be defined in the Heisenberg representation, 
i.e. satisfying the equation duk/dt = ih~^[HQ,ak]- To find the equation for the local Heisenberg 
operator 

am{l, t) = ^^y] ak{t) exp(iA;/) (2.10) 

V-/V , 

let us make the differentiation with respect to time Assuming the size of the space cell a is small, 
a <^ L, and a ^ 0, I may be considered as a continuous space coordinate. In this approximation 

d 

kex.p{ikl) = —i— ex.p{ikl) (2-11) 

as result it follows that the transfer equation for the local operator is 

-g-^+c-ja^{l,t) = (2.12) 

This equation is valid in a coarse space scale with a characteristic size a providing the fulfillment of 



(2.11), in this case the replacement dw -^ a5{l — V) is possible. 

For a many-particle problem ( p. 12 ) may be represented as set of coupled equations due to the 



interaction of local operators. The Hamiltonian Hq appears to be nondiagonal with respect to the 
indices //'. The equation of motion takes the form 

^a^(0 = ifi-^ [Ho; am{l)] = "^ E ^rn{l, l')am{n (2.13) 

I' 

the coupling constant 

links the oscillator in the point I with all neighbors. However, due to its 6 -shape, the coupling appears 
to be significant only for two adjacent oscillators. In other words, in the free space the interaction 
between local oscillators is described in a coarse scale by a derivative with respect to /: 

-Ql<^rn{l) = ^ E ^™(^' 0«m(0 (2.15) 



the equation ( |2.15| ) is obtained by differentiation of ( ^.lOf) with respect to / with taking into account 
( 2.11J ) Consider the commutator of local operators taking for different times. For the free evolution, 



when the dynamics of operators is determined by the Hamiltonian Hq 

[am{l,t);al^.{l,t + T)] = 5m„,' exp(icmr)- 5a{T) (2.16) 

the function (2.17) 

A^ ^ 1 has a sharp maximum for r ^ 0. Its weight is concentrated in the vicinity of the order of 
a/c = Av^, so it may be considered as a delta- function when At > a/c. Notice that (5a (0) = Aum- 
The occurence of the scale temporal delta-function is connected to the use of coarse space scale with 
a characteristic size a, where the time interval will be At > a/c. 

Let us introduce the interaction picture. Consider the slowly changing part or envelope in the 
local operator am{l)' 

am{i, t) = Am{l, t) exp{-iuJmt + iml) (2.18) 

where uj^ = cm, Then the unitary transformations in (|2.6| ) and ( |2.7] ) will take the form 
1 

7 

fc~m 

afc(t) = —j= ^ Am{l, t) exp{i(a;fc - w^)* - i{k - m)l} (2.20) 

the evolution of the operators a^ and Am{l,t) is defined only by the Hamiltonian of interaction V 

^^ak{t) = ih-'[V{t),ak{t)] (2.21) 

- + c-JAm{l,t) = ih~\V{t),Amil,t)] (2.22) 

3 Three - dimentional case. Quasioptical approximation 

Consider a normalized volume L^ with a number of cells A^ = NiN2N3,Ni = N2 = N3 = l/a, where 
for simplicity the cell was chosen to be of cubic shape. Then in the initial formulas the evident 
replacements will be: /c — > k, ?7i ^ m,l — > 1, x — > r. For the local operator in the interaction picture 



Am{l, t) = —= ^ ak{t) exp{-i(u;fc - u:m)t + i{k - m)l} (2.19) 

V-/V , 



the following expression should be written instead of (2.19): 
1 



A„,(l, t) = ^=J2 ak(t) exp(-i(u;k - Wm)t + (k - m)l) (3.23) 



k~m 



where the dispersion law has the form Wq = c-\/(q, q), q = k.m. Let us find the equation of motion 
for local operators. Similarly to the one-dimensional case let us differentiate ( p. 23] ) with respect to 
time. Assuming 1 to be a continuously changing variable, we will find instead of ( ^.11 ) 



d Id' d 
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(q,q)exp(zql) ^ "^5^ " ^ ^If + M ^"^^^^^^ ^'-^'^ 



When summing over k ~ m we replace k^ ^ m ^ niz This approximation corresponds to a quasiplane 
wave of frequency ujm = cm, propagating along the z-axis. In as much as a coarse spatial scale was 
introduced, for which vector 1 may be considered as a continuous function of the coordinates, let us 
make a replacement 1 — > r{x, y, z). As a result the equation of motion takes the form 

2 /a2 q2\\ 

+ — A^(r) = -in-'[V, A^{r)] (3.25) 




2u!m \ dx dy 



where V is the Hamiltonian of intraction. The equation following from ( 3.25| ) for the free field {V = 0) 



is well known in the classical theory. It describes the light propagation with diffraction taken into 
account in a quasioptical approximation. 

In the equation ( ^.25 ) the Hamiltonian V describing the electromagnetic field interactions with 



the medium should be expressed in terms of local operators. As an example consider the Hamiltonian 
of the light interaction with atoms in the dipole approximation, which is often a basis for a variety of 
problems. It has the form 




V{t) = —i V" \i T^^k exp(— icJkt + i'kvA)dA{t) + h.c. (3.26) 



The equation (3.26) is given in the interaction picture, where cLa is the operator of the dipole moment 



for an atom located at a point ta- Changing to local operators with the aid of ( 2.20 ), where the 
three-dimensionality should be taken into account, we will use the following approximations. Let the 
packet or the local mode interacts with the atom as a whole. It means that inside the band Af^ all 
frequencies ujk ~ uJm- Replace the atom position r^ by the position of the cell where this atom is 
located. Then the summation over the atoms may be divided into a sum over the cells 1 containing 
atoms and a sum over the atoms inside the cell. Suppose di is the operator of the atomic dipole 
moment in the cell 1, then the equation ( ^.26 ) will take the form Q 




^ (*) = -^Y.\l ^^^'"(l) exp(-iu;mt + im\)di + h.c. (3.27) 

Here a new normalization volume a^ appears, while the summation is performed only over the cells 
containing atoms. The obtained Hamiltonian describes the elementary interactions of the local field 
oscillators or photons in cells / with atoms located inside. As a result chaiging to local oscillators in 



the interaction Hamiltonian ( p. 261) reduces to an ordinary replacement: ak — > Ani(l), k -^ m, r — > 1. 
In a number of cases effective interaction operators obtained, for example, by unitary transfor- 
mations of the starting Hamiltonian (( p. 26 ), are used for the description of multi - photon processes. 



Given below are two effective interaction operators for two - photon and parametric interactions in 
the presentation of local operators: 

^2 = X! X! /l2 Ami (1) Am2 {^)Sl exp(i(cJ21 - ^mi " '^ma)*) + h.C. (3.28) 

1 1,2 

H^ = 5]5]Gi23(l)Ami(l)Am,(l)At„3(l)exp(-i(a;mi+Wm2-^m3)i) (3-29) 

1 1,2,3 

+ h.c. 



^details of these and following calculations are given in [tll[ 



where /12, G123 = 5123 exp(i(mi + m2 — ms)!) - are the couphng constants. The Hamiltonian H2 de- 
scribes two-photon interaction of the modes having frequencies in the region of two-photon resonance: 
a;2i ~ i^ki +'^k2- The operator Si = |2)/(1| corresponds to the transition of atoms located at the point 
I from the lower to the upper working level. The process of parametric interaction of three waves in 
a transparent medium conforms to the Hamiltonian H^^. 

4 Integral of motion and statistics of the light 

The integrals of motion may appear in the problems of light propagation in a medium, for which it is 
natural to use local description. The existence of the integrals enables to examine some peculiarities 
of the light statistics transformation without going into the solution of dynamic equations. Consider 
as an example two nonlinear processes: two-photon and parametric interactions, which are described 
by the effective Hamiltonians (3.281) and (3.2S|). Both processes are multimode, involving all pairs and 



triplets of modes having frequencies related as 

Wmi + '^ma = ^21 (4.30) 

Wmi + Wma = ^mg (4-31) 

Here after we restrict ourselves by the case of one-dimensional medium assuming the light propagating 
along the z-axis. 

The pairs of modes mi and m2 connected by the conditions ( f4.30| ) and (4.31) and (4.31) will be 



referred to as conjugated. For them the difference of the photon number operators 

I = ni{l) - n2{l) (4.32) 

where Uj = Al^.{l)Ajyi{l), j = 1,2 commutate with the Hamiltonian 

[^2,3,/] =0 (4.33) 

Consequently / and any function of the form /(I) are integrals of motion. For a problem with 
boundary conditions / will be considered as a continuous coordinate replacing / -^ z. Introducing the 
new variables and assuming the light propagating along the z-axis: t' = t — z/c, z' = z. Then by virtue 
of ( 4.33| ) the difference of the number of photons in conjugated modes is conserved for any point of 



the medium: 

I{z,t)=I{0,t-z/c} (4.34) 

With the use of the obtained operator integrals of motion it is possible to find the peculiarities of 
the transformation of the mutual correlation of conjugated modes for the propagation in a medium. 
Thus all statistical properties described by the correlation functions of differented intensity, e.g. of 
the form {{I{z, t)Y^I{{z, t + t))"^) = K^'P''^'{z, t) are conserved in virtue of ( [4.34 ). From the standpoint 



of observation of direct interest is the lowest order correlation function and its Fourier image 

POO 

K{z,n)= {I{z,t)I{z,T))e^p{inT)dT (4.35) 



It is possible to measure the function ( 4.35| ) in a scheme with two photodetectors (Fig.|^), where 
the fluctuation spectrum of difference photocurrent i^(0) is observed. For this scheme 

i\n)=r]{ni+n2)+v'^KNin) (4.36) 




Figure 1: 

Here the first independent of frequency term is the shot noise. Index A^ indicates normal ordering 
of the field operators to describe a detector that responds at the absorption of photon. The rate of 
registration r] = qAvm, where q is the quantum efficiency of the detectors, which for simplicity are 
assumed to be equal. The presence of the local mode Avm = c/a in the expression for the photocurrent 
spectrum follows from the formulas of detection given in the presentation of local operators. This may 
be elucidated as follows. A wide-band photodetector is needed to describ the light detection. Therefore 
to registrate the local mode in front of the detector we have to place an optical filter of a bandwidth 
Afm . Then Ai^rn that determines the rate at which the photons fall on the detector is the bandwidth 
of the scheme. 

Assuming that at detector the local operators commute as free field operators ( p.l6D let us to write 
the correlation function ( [4.35| ) in the normal and time-ordered form K[Q,) = i^7v(^) + ('^i +'^2)/^^'™- 
Then taking into account the motion integral one finds 



1 1 

Kn{z, Q) + (ni + 712) = Kn{0, Q) + (nio + n2o) 

AUm AVm 



(4.37) 



where (nio) and (^20) are the input photon number. As result the expression for the photocurrent or 
noise spectrum ( 4.36| ) takes the form 



i^{z, ^) = Aum{ni + n2)g(l - q) + Aum{niQ + n2o)qiq " 1) + ^^(0) ^) 



(4.38) 



where 



'(0, n) = Avmq{nio + n2o) + {AUraqfKN{Q, 0) 



(4.39) 



is the input noise. 

It is seen from (4.38) that for an ideal photodetector {q = 1) the input and output noise spectra 
of the light are equal. This means that the correlation of conjugated modes is conserved during 
propagation in the medium. Suppose that the input state of the modes is nonclassical, so it is a 
quantum correlation for which the shot noise is suppressed in bandwidth Avm by a factor of 1 — (7: 



'(0, Afl) = Avraq{niQ + n2o)(l 



(4.40) 



The state of the light resulting in (4.40) produces, for example, by optical parametric amplifier (OPA). 
Then the level of the shot noise suppression in the output is unchanged: 



\z,AVt) = Av^q{ni+n2){l 



(4.41) 



This means that the quantum property of the hght in the medium is conserved. Moreover in the 
medium with two-photon or parametric interaction the conjugated modes may be amphfied, at that 
the observed initial correlation, in particular the quantum correlation, is conserved. This specific 
feature determines the properties of spontaneous radiation. Thus, if input is the vacuum state, then 
i'^{z, AQ) = at the output. This means that the conjugated modes arising in spontaneous radiation 
have nonclassical correlation, that results in a suppression of shot noise. Indeed all these properties 
are following only from the existence of integrals of motion in systems expanded in space. 

5 Local quasiprobabilities 

In previous items the Heisenberg picture was introduced for local operators. In the Schrodinger picture 
the evolution is determined by the density matrix of electromagnetic field p, which may be connected 
to c-number functions P{{ak}, s) called s-ordered quasiprobabilities. They arise in the density matrix 
expansion over operators A({afc},s), which are Fourier-images of s-ordered displacement operators 
forming the complet set: 

A({a,.},s) = - /'{d2/3jj]exp{s|/3fc|2 + (/?*(«;, -a^)- /i.e.)} (5.42) 

^ -^ k 

{cPPk} = Hkd'^Pk- Then the s-ordered quasiprobability is determined by the expression 

P{{ak};s) = Sp{A{{ak},s)p) (5.43) 

Using the replacements a^ — > am{l) in the above formulas we obtain the local quasiprobability 
-P({am(0}i s). In the following we will restrict ourselves by the case s=l corresponding to the nor- 
mal ordering of the field operators which is described by the Glauber-Sudarshan quasiprobability 
P{s = 1) = P. This function arises in the density matrix expansion over coherent states or diagonal 
representation. Prom the given expressions it is possible to obtain the relation between local and non- 
local quasiprobabilities. Similary the case of operators, it is possible to proceed from one distribution 
function to another by the change of variables: 

P{{o^mil)}) ^^ nWk}) (5.44) 

{amim ^^ {ak} (5.45) 



where the variables are related by the expression of the (2^) and (|27 

In the Schrodinger representation the peculiarities of the local description specified by the many- 
particle character of the problem become apparent to the same extent as in the Heisenberg represen- 
tation. Here the transfer equations occur, which are formally equivalent to the BBGKI chains for the 
partial distribution functions. Of all the hierarchy of the distribution functions in the present case 
the following two types of one-particle quasiprobabilities are of interest. They arise from the function 
-P({am(0})i of all modes of all local oscillators. Averaging over all oscillators except for the chosen 
one we will find the distribution function P{af^-\{1)) for all modes of one local oscillator located at 
a point /. The averaging of P(ar„i (Z)) over all modes except for one results in a function P{am{l)) 
describing one mode of one local oscillator. 

For a free evolution the introduced one-particle distribution functions Pi = Pia{m}{^))^ -P(cim(0) 
obey the following transfer equation: 



In ( ^.461) the derivative over / describes the interaction of local oscillators or excitation transfer. 
Similarly to ( |2.13| ) it may be presented in the form 

Sp{a^{l)) = -hm{l, l)P{am{l)) - E / d^oim{l')hm{h l')P{am{l), a^{l') (5.47) 

where hmil-,1') = i^mil-,1 ){d / dam{l))(y.m{l ) + c.c. is the differential operator of free evolution with 
Hamiltonian Hq <-> Yl,mU' hm{l, I')- 

6 Fokker - Planck and Langevin equations 

An approach based on the master equation for the quasiprobability P{{ak}) in the Fokker-Planck 
approximation is often used to describe the statistical properties of light. In the problems of the 
light interaction with atoms the master equation for the electromagnetic field can be obtained by an 
adiabatic elimination of atomic variables. It is not necessary to derive the field equation once again 



in term of the local description, because one may immediately use the change of variables (5.45). 

As an example consider the interaction of light with a two-level system in the lowest approximation 
for which the processes of the linear amplification or absorption type occur. Such a medium is described 
by a linear susceptibility, K{io) = \d\'^{uJo — lo — ij)^^ where d and luq is the dipole moment and the 
frequency of the atomic transition, 7 is the decay rate or transverse relaxation The form of the field 
equation obtained by the adiabatic elumination of the fast atomic variables is well known (see, 0): 



d_ 
dt 



PiM) = ^ (^w^». + '^'"^iSh + "-') ''""'» *"'* 



Here the coefficients are defined by the linear susceptibility and populations of the upper A'^2 and lower 
A^i levels: 

Aik) = CkiNi - N2)Im K{ujk) - iek{Ni + N2)Re ^(wfc) 
Q{k) = ekN2lm K^uJk) 
€k = h-'^{hu;k/2£oL^) 

When changing the variables we will assume that the wave packet forming the local mode interacts as 
a single whole. Then inside the band A; ~ ttt. it is possible to neglect the dispersion of all harmonics. 
This means that A(k) ~ A{m), Q{k) ~ Q{m). This approximation enables us to write at once the 
transfer equation for one-particle quasiprobability, say P(a{m}(/)) = Pi in the form 



^ " ^^) ""' = ? (■'""'*£(o"""> ^ «'">ra<(o ) ^' ^ "■ "'■^'' 



Notice that the structure of differential operators in right- hand site is the same in ( 6.49 ) and ( 6.4S| ), 
and the whole of the passage procedure to the local description reduces to an addition of a derivative 
over / to the left-hand side of the equation. 

The equation of quasiprobability ( |6.48| ) is in agreement whith the following Langevin equations: 

^afc = -A{k)ak + fk (6.50) 



d_ 
dt 



al = -A*{k)al + fl 
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where the correlator of random forces is defined by a diffusion coefficient 

if kit) flit + r)) = 2Q{k)6kk'S{T) (6.51) 

In as much as the Langevin variables here are associated with the diagonal representation they corre- 
spond to the normally ordered averages of the field operators {a\{t)ak{t + r)) = {a%{t)ak{t + r)). 

In the local description the main moment in the formulation of Langevin equations is the finding 
of the random force correlator. To this end we use changing of variables of the form ( |2.19| ), (|2.20| ), 
where ^m(0 ~^ CKm(0 ^k ~^ Cik-, thereby we introduce the local Langevin variable am(0- Then the ap- 
proximation A[k) K, A{m)^ Q{k) ~ Q{m) used in writing down the equation for local quasiprobability 
leads to Langevin equations 

^ + c^jam(/,t) = -A{m)am{l,t) + fm{l) (6.52) 

d d 



(^ + c^)<G,t) = -A*im)a*^il,t) + flil) 



For the used variables there follows an expression for a random source 
1 

7 

fc~m 



fmil, t) = -j= V fk ex.p{-i{ujk - iOm)t + i{k - m)l) (6.53) 

V-'V I, 



From this we will find that the correlator of random sources as well as in ( |6.51| ) is defined by the 
diffusion coefficient 

{f^il\t)fl{l,t + T)) = 2Q{m)6^^^5n^5{T) (6.54) 

As to ( |6.54D it should be noted that during the transition to a local description with a coarse space scale 
a scaling time (5- function should be used defined as in ( 2.17] ) and the / magnitude should be considered 



as a continuously coordinate. This implies that in ( |6.54| ) the replacements take place 5{t) — > ^^(t) 
5ii' -^ a6{l — V). As a result the correlator of random sources in local Langevin equations has the form 

{fm'iit)fUt + r)) = 2Q{m)6mm'a6{l - r)<5„(r) (6.55) 

In comparison with the ordinary description in ( |6.55D there arises a space delta-function due to the 
local correlation of oscillators in space. 

7 Diffraction in linear medium 



The differential operator in ( ^.25 ) becomes non-Hermitian due to diffraction effects, so the correspond- 



ing equation for quasiprobability becomes very complicated. However it is possible to use the Langevin 
formulation. In this way for the case of a linear medium considered above the change of variables in 
nonlocal Langevin equations ( |6.52|) leads to Langevin equations which in the three-dimensional case 
have the form 

a~ + a~ "»" = -A{m)am + fm{r,t) (7.56) 




2u!m \dx dy 
l + "^+^2^(l^ + |))"^ = -A*{m)a*^ + fL{r,t) 
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where the random force correlator is defined as 

(/m'(r',t)/^(r,t + r)) = 2Q{m)6mn.'a^6{r - r')6a{r) (7.57) 

Here Q{m) is the diffusion coefficient, which may be obtained from the corresponding ( [7.56D but 
nonlocal Fokker-Plank equation ( |6.48| ): Q{m) ~ Q{k). 

The equations (7.56) are linear in field amplitudes, so they are easily integrated. Let us denote the 



transverse radius-vector s = (x, y) and introduce a travelling coordinate system t' = t — z/c, dz' = z. 
Then for the given conditions at the boundary one finds 

A 



A f 
am{z,s,t) = exp( z) (fsi am{0,si,t — z/c) U{Osi\zs) +W{z,s,t — z/c) (7.58) 

W = - dzi exp( (z — zi)) cfsifm(zi,si,t — z/c)U(ziSi\zs) (7.59) 

c Jo c J 



- / azi exp( 

c Jo c 

where the Green function is defined as 



U(ziSi\zs) = —i — rexp im (7.60) 

27r(z — zi) \ z — zi I 

Here am(0, si, t — z/c) is the transverse field distribution in the input, A = A{in) 
Consider the properties of a random source W, defining its correlation function 

{Wiz, 0, t)W\z, s, t + T)) = ^Qa^Sair) T exp f-^-^{z - zi)) D{z - zi)dzi (7.61) 



where the value 

m 



D 



2tt{z — z\) 




Tfl 

d^si exp ( i sis 1 (7.62) 

z — Z\ 



describes the diffraction effects. Since the diffusion coefficient Q is determined by the number of atoms 
at the upper level A'^2 the transverse sizes of the medium where the light propagates should be taken 
into account. In this way the limits of integration over si in D are defined. Thus D = 5{s) if the 
medium is not limited in the transverse direction. This means that the noise source is (5-correlated and 
the noise is white. Taking into consideration the finiteness of the size in transverse direction changes 
the situation. Thus in a medium confined by a cylinder with a radius R 

mR ( msR \ ( . ms^ \ 

where J\ is the Bessel function. 
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